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Abstract 

A comprehensive universal description of the rotational-vibrational spectrum for two identical 
particles of mass m and the third particle of mass m\ in the zero-range limit of the interaction 
between different particles is given for arbitrary values of the mass ratio m/m\ and the total 
angular momentum L. It is found that the number of vibrational states is determined by the 
functions L c (m/mi) and Lb(m/mi). Explicitly, if the two-body scattering length is positive, the 
number of states is finite for L c (m/m\) < L < L\,{m/mi), zero for L > L^m/mi), and infinite 
for L < L c {m/m\). If the two-body scattering length is negative, the number of states is zero 
for L > L c (m/m\) and infinite for L < L c (m/m\). For the finite number of vibrational states, 

all the binding energies are described by the universal function s ln (m / mi) = £(£,77), where 

N — 1/2 / m 
£ = — =, r? = a ; r, and N is the vibrational quantum number. This scaling 

y^(L + l) V m i L ( L + 1 ) 

dependence is in agreement with the numerical calculations for L > 2 and only slightly deviates 

from those for L = 1,2. The universal description implies that the critical values L c (mjm\) 



and Lj,(m/mi) increase as 0.401ym/mi and 0.563 y / rn7mi, respectively, while the number of 
vibrational states for L > L c {mjm\) is within the range N < N max ~ 1.1-v/ L(L + 1) + 1/2. 



PACS numbers: 03.65. Ge, 21.45. +v, 03.75.Ss, 36.90.+f 
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I. INTRODUCTION 



The universal low-energy few-body dynamics of two-species compounds is of much interest 
both for atomic and many-body physics. In this respect, the study of the three-body energy 
spectrum gives insight into the role of triatomic molecules and few-body scattering. The 
area of applications includes the investigation of multi-component ultra-cold quantum gases, 
e. g„ binary Fenni-Bose and Fermi fl Q, mixtures and of impurities embedded in a 
quantum gas Ub|. which are present.y under thorough experiments and theoretic*, study. 
In addition, one should mention the reactions with negative atomic and molecular ions [8, 9J. 
The universal isotopic dependence of the three-body energy spectrum was multiply dis- 



cussed 
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14|, nevertheless, the main objective was the description of Efimov's 



spectrum. Recently, the infinite number of the 1 + bound states was predicted [15( for three 
identical fermions with the resonant p-wave interaction. Concerning the low-energy scatter- 
ing, one should mention a two-hump structure in the isotopic dependence of the three-body 



181 ] and the two-component model 



recombination rate of two-component fermions 
for the three-body recombination near the Feshbach resonance 

The main aim of the paper is a comprehensive description of the finite three-body 
rotational-vibrational spectrum in the zero-range limit of the interaction between differ- 
ent particles. Both qualitative a nd numerical results are obtained by using the solution of 
hyper-radial equations (HREs) 120. 21 



221 ]. The detailed study of the bound states and 



scattering problems for the total angular momentum L = 1 was presented in [18]. 



II. OUTLINE OF THE APPROACH 



Particle 1 of mass mi and two identical particles 2 and 3 of mass m are described by 
using the scaled Jacobi variables x = y/2jl (r2 — ri) , y = y/2p[r^ — (miV\ + mr2) / (mi + m)\ 
and the corresponding hyper-spherical variables x = pcosa, y = psina, x = x/x, and 
y = y/y, where is the position vector of the ith particle and \x = mmijim + mi) 
and ft = mim + mi) /(mi + 2m) are the reduced masses. In the universal low-energy 
limit, only the s-wave interaction between different particles will be taken into account 
provided the s-wave interaction is forbidden between two identical fermions and is strongly 
suppressed between two heavy bosons in the states of L > 0. The two-body interaction 



2 



is defined by imposing the boundary condition at the zero inter-particle distance, which 
depends on a single parameter, e. g., the two-body scattering length a [3]. This type of 



interaction is known in the literature as the zero-range potential 



the Fermi 



Fermi-Huang 



2J] or 



251 ] pseudo-potential, and an equivalent approach is used in the momentum- 



space representation 



261 ] . The units H = 2// 



1 are used throughout; thus, the 



binding energy becomes the universal function depending on the mass ratio m/rax and the 
rotational-vibrational quantum numbers L and N. In view of the wave-function symmetry 
under permutation of identical particles, a sum of two interactions between different particles 
is expressed by a single boundary condition at the zero distance between particles 1 and 2, 



lim 



d a 

— tana — p— 

oa a 



* = . 



(1) 



The problem under study is conveniently treated by using the expansion of the properly 
symmetrized wave function, 



[1 + SP) h ^ y l f>(p)^(a,p) 

n=l 



p5/ 2 sin 2a 



(2) 



which leads to the hyper-radial equations for the functions f n (p) |18(]. Here P denotes 
permutation of the identical particles 2 and 3, S = 1 and S = — 1 if these particles are 
bosons and fermions, respectively, YlmIy) is the spherical function. The action of P on the 



angular variables in the limit a — > ir/2 is given by PYlm{y) 



-l) L Y LM (y) and Pa 



where uj = arcsin(l + mi/m) 1 . The functions <£>^(a, p) in the expansion ([2]) are the solutions 
of the equation on a hypersphere (at fixed p), 

d 2 L(L + 1 



da 2 



sin 2 a 



+ 7n(p) 



<P n (<*,P) = 



(3) 



complemented by the boundary conditions ip^(0,p) = and 



a-n/2 \aa \a\ J sm2o; 



(4) 



where a set of discrete eigenvalues l 2 {p) plays the role of the effective channel potentials in 
a system of the hyper-radial equations 18|] . The functions satisfying Eq. ([3]) and the zero 
boundary condition are straightforwardly expressed 27| via the Legendre function 



(5) 
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The functions 0x )7 (a) are odd functions on both variables 7 and a satisfying the recurrent 
relations sin a 0£ + i i7 (a) = (7 — L — 1) cos a 0L 7 (a) — (7 + L)</>£ i7 _i(a), which follow from 
those for the Legendre functions. It is convenient to write i;7 (a) = A i;7 (cot a) sin 7a + 
i?L j7 (cota) cos 7a, where Al^(x) and B^^ix) are simple polynomials on 7 and x, which are 
explicitly given for few lowest L by A 0)7 (x) = 1, -B 0)7 (x) = 0, A ln (x) = —x, Bi r/ (x) = 7, 
A 2}1 {x) = 1 - 7 2 + 3x 2 , B 2a (x) = -37X, A 3)7 (x) = 3x(27 2 - 3 - 5x 2 ), and S 3i7 (x) = 
7(15x 2 + 4 — 7 2 ). Substituting ([5]) into the boundary condition (TJJ and using the identity 
0L+i, 7 ( 7r /2) = ^da^ one comes to the transcendental equation for 7 2 (p), 

a=n/2 

p± L7 (tt/2) = i+ll7 (7r/2) - ^zl!0 Li7 H . (6) 
|a| sin 2a; 

The attractive lowest effective potential determined by 7 2 (p) plays the dominant role 
for the binding-energy and low-energy-scattering calculations, while the effective potentials 
in the upper channels for n > 2 contain the repulsive term 7^(p)/p 2 and are of minor 
importance. Thus, a fairly good description will be obtained by using the one-channel 
approximation for the total wave function (T5]) where the first-channel radial function satisfies 
the equation 



r d 2 7i 2 (p) - 1/4 



+ E 



fi(p)=0. (7) 



dp 2 p 2 

Note that the diagonal coupling term is omitted in Eq. (CO), which does not affect the final 
conclusions and leads to the calculation of a lower bound for the exact three-body energy. 
Our calculations [3] shows that the one-channel approximation provides better than few 
percent overall accuracy of the bindin g en ergy. 



The most discussed feature 



11 



12 



13| of the three-body system under consideration 



is the infinite number of the bound states for small L and large m/m\ (more precisely, for 
the finite interaction radius the number of states unrestrictedly increases with increasing 
\a\/r ). As the effective potential in <£7§ is approximately given by (7 2 (0) — 1/4) / p 2 at small 
p, the number of vibrational states is finite (infinite) if 7 2 (0) > (7 2 (0) < 0). According 
to Eq. (J6]), 7 2 (0) decreases with increasing mjm\ and becomes zero at the critical value 
(m/mi)^. Thus, one can define the step-like function L c (m/mi), which increases by unity 
at the points (m/rrii)^, so that the number of vibrational states is infinite for L < L c {mjm^) 
and finite for L > L c (mjm\). Solving Eq. (jEJ) at 71 — ► and p — > 0, one obtains the 
exact values {m/mi) cL) which approximately equal 13.6069657, 38.6301583, 75.9944943, 
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125.764635, and 187.958355 for L = 1 — 5. Originally, the dependence L c {mjm\) was 
discussed in |lo| . 

Analyzing the eigenvalue equation ([6]) one concludes that for a > and S(— ) L = — 1 
the effective potential exceeds the threshold energy E = —1, 7f(p)/p 2 > —1, therefore, the 
bound states only exist if either two identical particles are bosons and L is even or two 
identical particles are fermions and L is odd. Furthermore, one obtains the trivial answer if 
a < and L > L c (m/mi), for which 7^(p) > and there are no three-body bound states. 



III. NUMERICAL RESULTS 



The mass-ratio dependence of the binding energies ££jv(m/mi) for L > L c {mjm\) and 
a > is determined numerically by seeking the square- integrable solutions to Eq. (j7j). 
Mostly, the properties of the energy spectrum are similar to those for L = 1, which were 
carefully discussed in 18|]. For given L, there is the critical value of mjm\ at which the first 
bound state arise, in other words, there are no three-body bound states for L > Lb(m/mi), 
where the step-like function L h {m/mi) undergoes unity jumps at those critical values. Fur- 
thermore, all the bound states arise at some values of m/mi being the narrow resonances 
just below them. For the mass ratio near these values, the binding energies and resonance 
positions depend linearly and the resonance widths depend quadratically on the mass-ratio 
excess. Exactly at these values one obtains the threshold bound states, whose wave func- 
tions are square-integrable with a power fall-off at large distances. A set of these values of 
m/mi (more precisely, the lower bounds for them) is obtained numerically and presented in 
Table [H With increasing m/mj, the binding energies monotonically increase reaching the fi- 
nite values (shown in Tabled]) at (m/m 1 ) cL ; just below (m/mi) c i they follow the square-root 
dependence on the difference m/rrii — (mjm^)^. Correspondingly, the number of the vi- 
brational states increases with increasing mjm\ taking the finite number N max at (m/mi) c i 
and jumping to infinity beyond (m/m 1 ) c i; in the present calculations N max — L + 1 for 
L < 9 and N max = L + 2 for 10 < L < 12. 
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TABLE I: Upper part: Mass ratios for which the Nth. bound state of the total angular momentum 
L arises. Lower part: Binding energies eln for the mass ratio fixed at (m/mi) C L- 

N L = 1 L = 2 L = 3 L = 4 L = 5 



1 7.9300 22.342 42.981 69.885 103.06 

2 12.789 31.285 55.766 86.420 123.31 

3 - 37.657 67.012 101.92 142.82 

4 - - 74.670 115.08 160.64 

5 - - - 123.94 175.48 

6 - - - - 185.51 

1 5.906 12.68 22.59 35.59 52.16 

2 1.147 1.850 2.942 4.392 6.216 

3 - 1.076 1.417 1.920 2.566 

4 - - 1.057 1.273 1.584 

5 - - - 1.049 1.206 

6 - - - - 1.045 



IV. UNIVERSAL DESCRIPTION OF THE SPECTRUM 

A comprehensive description of the spectrum is obtained by using the large- L (correspond- 
ingly, large-m/mi) asymptotic expression for the binding energies £LAr(m/mi). Taking the 
quasi-classical solution of Eq. (Ej) satisfying the zero boundary condition, 

_ — \ L/2+1/4 

Vl + x 2 sin a — cos a \ , . 

; = ,j = , (8) 

V 1 + x sin a + cos a J 

where 71 = in and x = k/(L + 1/2), one writes the eigenvalue equation (jBj) in the form, 



P 



2 exp (k arcsm j / Vl + x 2 sm 2 u - cos u \ 



\ L/2+1/4 



L + l/2 (L + l/2)sin2^ \Vl + x 2 sin 2 uj + cos u J KJ 

In the limit of large L and mjm\ the eigenvalue equation (jUj) reduces to 

pcosu; = w — e~" , (10) 

where u = cosu y/n 2 + (L + 1/2) 2 . Notice that taking the limit k -> and p -> in ([TO]) 
one immediately obtains the relation coscij cL = u /(L + 1/2), where sino; ci = (m/m 1 ) ci /[l + 



(m/mi) c i] and Uq ~ 0.567143 is the root of the equation result, one finds the 

asymptotic dependence {m/mi) c L ~ 6.2179(L + 1/2) 2 and the inverse relation L c {mjm\) w 
{uQyj2m/m\ — l)/2 « 0.40103A/m/mi — 1/2. Now the asymptotic dependence £LN(m/mi) 
for large L and m/mi can be obtained by the quasi-classical solution of ((71) with jf(p) = 
(L + 1/2) 2 — [u(p)/ costu] 2 and u(p) determined by (TTUl) . 

■+ i_l e -u 

du Jv? - e LN (u - e- u ) 2 - L(L + 1) cos 2 u = ir(N - 1/2) cosu , (11) 

u — e u 

where U- and u + are zeros of the integrand. 

Following Eq. (TTTj) . one expects to express the binding energies via the universal function 

£LN( m / m i) = °f two scaled variables £ = — = = and n = ,/ This 

a/£(£ + 1) ymiL(L + l) 

two-parameter dependence is confirmed by the numerical calculations (up to m/mi ~ 700), 
which reveal that the calculated energies for L > 2 lie on a smooth surface as shown in 
Fig. [TJ Even for the smallest L — 1,2 the calculated energies are in good agreement with 
the two-parameter dependence showing only a slight deviation from the surface. 



(-E) 



1/2 



0.75 



0.5 



0.25 




FIG. 1: Universal dependence of the bound-state energy E on the scaled variables £ 



N — 1/2 
VHL + 1) 



and ?7 



m 



-. The calculated values for L = 3, 12 are plotted by symbols. The surface 



m 1 L(L + 1) 

boundary and its projection on the £ — 77 plane are shown by solid lines. 
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The variables £ and 77 take values within the area limited by the line £ = 0, the line 

Tj = T] max « a/2/mo ~ 2.493574 stemming from the condition of finiteness of bound states 

L > L c (m/mi), and the line £(£,77) = 1 expressing the condition of arising of the bound 

states at the two-body threshold. As shown in Fig. [T], the smallest value 77 = f]mi n is 

at £ = 0, which corresponds to the condition of arising of the first bound state in the 

large-L limit. To find it one requires that u + = w_ = Ub at eln = 1 in Eq. (fit)]) , which 

leads to r] min m a/2/(m^ — 1) ~ 1.775452, where Ub ~ 1.278465 is the root of the equation 

u = 1 + e~ u . This gives the asymptotic dependence for arising of the first bound state, 

Lb ~ V m lnV m / mi — 1/2 ~ 0.563237A/m/rrii — 1/2. At the line 77 = r] ma:r the variable £ 

takes its largest value £, max , which determines the large-L dependence of the number of the 

vibrational states N max for a given L. The calculation of the quasi-classical integral ffTTj) 

gives m_ = w , u + « 2.872849, £ maa; = / -v/e u (2?z - w^m) _ 1 du w 1.099839, 

7TWo A. ^e" — 1 v 

and the large-L estimate A^ max = £ max ^L(L + 1) + 1/2. Taking the entire part of this 
expression, one can predict that the dependence N max = L + 1 for L < 10 changes to 
N m ax = L + 2 at L = 10, which is in agreement with the numerical result. 

The universal surface £(£,77) is bound by three lines, which are described by fitting the 
calculated energies for L > 3 to simple dependencies plotted in Fig. HI As a result, the 
line defined by £(£,77) = 1 is fairly well fitted to r] = (i] min + a£)[l — c£(£ — £ max )], where 
a = {Vmax-Vmin) / t,max ~ 0.652933 is fixed by the evident condition S(C, max , f]max) = 1 and the 
only fitted parameter is c = 0.1865. Furthermore, the analysis shows that the next boundary 
line defined by r] = r] max is described by £~ 1/2 (£,r] max ) = a^(l - a 2 £)[l - c i£(£ ~ Cmax)}, 
where a% — 1 + Uo ~ 1.56714 is fixed by the asymptotic behaviour of the integral (fTTj) at 
L — > oo and 77 = r] max , a 2 ~ 0.38171 is fixed by the condition ai^ max {l — a 2 £ ma z) = 1, and 
the only fitted parameter is C\ = 0.1881. In particular, at the critical mass ratio the binding 
energy of the deep states in the limit of large L is described by £((,,r] max ) — > (ai£)~ 2 , i. e., 
£jVi[(W m i)ci] — (N-i/2)^(i > +uo) i ■ T ne ^ird boundary line at £ — >• is described by the 
dependence £~ 1/2 (0,r]) = a 3y /r] max - rj[l + c 2 (i] - r) min )), where a 3 = 1/ y/rjmax ~ Vmin ~ 1-18 
is fixed by £(0,f] m i n ) = 1 and the only fitted parameter is c 2 = 0.3992. 
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V. CONCLUSION 

The presented results complemented by the accurate calculations for L = 1 [18] provide 
in the universal low-energy limit a comprehensive description of the rotational- vibrational 
spectrum of three two-species particles with the short-range interactions. Essentially, all 
the binding energies are described by means of the universal function £(£,77) for those 
L c (m/mi) < L < L^imjm^) which correspond to the finite number of vibrational states. 
One expects that the universal picture should be observed in the limit \a\ — > oo, e. g., i f the 
potential is tuned to produce the loosely bound two-body state as discussed in |22|, |28 |. 

It is of interest to discuss briefly the effect of the finite, though small enough interaction 
radius r <C a- For L < L c (m/mi) Efimov's infinite energy spectrum is extremely sensitive 
to the interaction radius tq and to the interaction in the vicinity of the triple-collision 
point, whereas for L > L c (m/mi) the binding energies depend smoothly on the interaction 
parameters provided r <C a. For this reason, one expects not an abrupt transition from the 
finite to infinite number of bound states for L = L c {mjm\) but a smeared off dependence 
for any finite value of r /a. 

It is worthwhile to mention that arising of the three-body bound states with increasing 
mass ratio is intrinsically connected with the oscillating behaviour of the 2 + 1 elastic- 
scattering cross section and the three-body recombination rate. In particular, for L — 1 it 
was shown in 18j that two interference maxima of the scattering amplitudes are related to 
the arising of two three-body bound states. Analogously, the dependence of the scattering 
amplitudes on the mass ratio for higher L would exhibit the number of interference maxima 



which are related to arising of up to N max = 1.099839 V /L(L + 1) + 1/2 bound states. 

Concerning possible observations of the molecules containing two heavy and one light 
particles in the higher rotational states, one should mention the ultra-cold mixtures of 87 Sr 

n 

with lithium isotopes [18( and mixtures of cesium with either lithium or helium. In particular, 
for 133 Cs and 6 Li the mass ratio m/mi « 22.17 is just below the value m/mi = 22.34 at 
which the L = 2 bound state arises and m/rrii « 33.25 for 133 Cs and 4 He is above the value 
m/rrii = 31.29, which corresponds to arising of the second L = 2 bound state. Also, a 
complicated rotational- vibrational spectrum and significant interference effects are expected 
for the negatively charged atomic and molecular ions for which the typical total angular 
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momentum up to L ~ 100 becomes important due to the large mass ratio. 
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